Thick films coating a plate withdrawn from a bath 



J.H. Snoeijer 1 , J. Ziegler 1 , B. Andreotti 2 , M. Fermigier 2 , J. Eggers 1 
School of Mathematics, University of Bristol, University Walk, Bristol BS8 1TW, UK 1 
PMMH, UMR 7636 CNRS-ESPCI-P6-P7, 10 rue Vauquelm, 75231 Pans Cedex 05, France 2 

(Dated: March 27, 2008) 

We consider the deposition of a film of viscous liquid on a flat plate being withdrawn from a 
bath, experimentally and theoretically. For any plate speed U, there is a range of "thick" film 
solutions whose thickness scales like U 1 ^ 2 for small U . These solutions are realized for a partially 
wetting liquid, while for a perfectly wetting liquid the classical Landau-Levich-Derjaguin (LLD) film 
is observed, whose thickness scales like U 2 ^ 3 . The thick film is distinguished from the LLD film by 
a dip in its spatial profile at the transition to the bath. We calculate the phase diagram for the 
existence of stationary film solutions as well as the film profiles, and find excellent agreement with 
experiment. 



When a solid object is pulled out of a liquid reser- 
voir, a thin layer of liquid is entrained by viscous drag 
(Fig. [TJ . This principle is used widely in coating technol- 
ogy, where it is known as dip coating, because it is one of 
the simplest ways to deposit a thin film of liquid on a sub- 
strate pp. According to the pioneering work of Landau 
& Levich [5] and Derjaguin [3] (LLD), a film of unique 
thickness Iilld is selected by the speed of withdrawal U. 
The LLD solution has remained the basis of coating the- 
ory for more than 60 years, having been generalized to 
include the effects of inertia |4j |5] , deposition on curved 
substrates [6] , and non-Newtonian fluids [7] . The relative 
size of viscous drag and capillary retention in the film is 
measured by the capillary number Ca = Utj/j, where i] 
is the viscosity and 7 the surface tension. At the foot of 
the film, LLD introduced the requirement that the film 
be matched smoothly to a static capillary meniscus (see 
Fig. [TJ, whose size is controlled by the capillary length 
= \J II 1 P9- In the limit of small Ca, this matching 
yields Iilld = 0.9464Ca 2/3 , which gives the small-Ca 
behavior of the LLD line in the phase diagram, cf. Fig. [2] 

In this Letter, we show that at a given speed there 
exists another, thicker film in addition to the LLD film 
(cf. Figs Jl] and |2). These new solutions do not match 
smoothly to the bath, but exhibit a bump at the foot of 
the film (cf. Fig. [l](c)). As we will explain in more detail 
below, experimentally the thick solutions are most easily 
produced for partially wetting fluids just above the criti- 
cal plate speed Ca/ at which a film begins to be deposited 
[5]. Our theoretical analysis however is for the long-time 
limit only, in which the plate is covered completely by a 
film, so any solution is characterized by a pair (/ij,Ca), 
as shown in Fig (2] 

A striking feature of the new solutions is that they 
cover a continuous range of thickness hf (Fig. [2] grey 
area) . Each of these solutions entrain a different amount 
of liquid flux q, representing the volume of liquid with- 
drawn from the bath per unit time and unit plate length. 
Clearly, this flux is limited by the flow through the nar- 
row dimple. Below we compute the maximum possible 
flux through the dimple and show how this provides a 
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FIG. 1: Fluid films being deposited on a plate withdrawn 
from a bath of viscous liquid; symbols are experimental 
measurements of the film profile at successive times, (a) 
Schematic of the experiment, (b) The fluid wets the plate 
and deposits a LLD film, Ca = 9.27 x 10 -3 ; the solid line is 
the classical LLD solution, which predicts hf = 57pm. (c) 
The plate is treated such that the fluid wets partially, and 
deposits a thick film, Ca = 8.05 x 10 -3 ; the solid line is the 
result of our theory with the film thickness fitted to the ex- 
perimental result of hf = 198rtm. 

lower bound on hf (Fig. |2j. For larger hf, draining due 
to gravity becomes increasingly important, so at the up- 
per bound (Fig. [2] dashed line) the flux vanishes. Above 
the dashed line, no solutions are possible, unless the film 
is alimented from above. 

In our theoretical analysis we rescale all lengths by the 
capillary length £ c . The interface profile is analyzed using 
the lubrication equation [5] for free surface flows [TO] : 

«'-l + ^(Ca-i)=0. (1) 
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FIG. 2: Phase diagram of stationary solutions in the parame- 
ter space (hf, Ca). The dimple solutions exist in the light grey 
region, but only those in the dark region actually entrain liq- 
uid from the bath. The dotted line shows the Landau-Levich 
solution. The open symbols are for a partial wetting situ- 
ation, the full symbols for complete wetting; error bars are 
below the size of the symbols. The two round symbols corre- 
spond to the two measurements shown in Fig. [l] The arrows 
show the experimental path as the speed is slowly increased 
starting from the dashed line. 



Here h(z) is the interface thickness, k the interface cur- 
vature and q the flow rate. The three terms correspond 
to the capillary pressure, gravity, and viscous effects, re- 
spectively. For steady solutions, considered here, q is 
constant along the plate. Film solutions become asymp- 
totically flat, (k = 0), so that the flux can be related to 
the film thickness hf as 



that the curvature h" remains finite in the limit of small 
Ca, which we are analyzing. This ensures that (|3| can 
be matched to the capillary meniscus of the bath, which 
has a constant curvature of \f2 so the boundary 

condition in the limit rj — ► — oo is H"^ = \/2- Defining 
Q = gCa -5 / 3 , one obtains the equation for the dimple 
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(4) 



So far Q is the same as for the LLD solution, but the 
boundary condition toward the flat film region will be less 
restrictive, so we are able to produce additional solutions. 
In that region, the curvature becomes small, giving the 
condition H'^ = 0. Below we will justify this argument 
in more detail. 

Figure [4] shows that there is a one-parameter family of 
solutions H(rj), each corresponding to a different value 
of Q, and forming a dimple that matches the thick film 
to the static bath. The case Q — can be solved analyt- 
ically [12], while the others have been obtained numer- 
ically. We identified a maximum value Q ma x ~ 1-376, 
above which the boundary conditions can no longer be 
satisfied. Note that the LLD film is a particular solution 
of Q, for which H approaches a constant value rather 
than shooting back upwards (dotted line). 

The final step is to relate the flux through the dimple to 
the film thickness hf at a large distance from the reser- 
voir. Just before q{hf) changes sign, it becomes small 
enough to match the flux in the dimple (cf. Fig. |3j. For 
a small enough capillary number, ^ can be solved to 
give 



h f ~ (3Ca) 



1/2 



-QCa ; 



2/3 



(5) 



q= hf Ca 



7 



(2) 



Crucially, q exhibits a non-monotonic dependence on 
h f , owing to gravitational draining (Fig. [3Jd) . For small 
h f the fluid velocity is uniform across the layer and equals 
the plate velocity, yielding a (dimensionless) flux /i/Ca. 
For larger thickness gravity induces draining inside the 
film, which can even reverse the direction of liquid trans- 
port: q < corresponds to transport from the film into 
the reservoir. As a consequence, the same flow rate q 
can be supported by two very different thicknesses of the 
film. This is why, physically, the thin 'dimple' region can 
be matched to the much thicker film (Fig. [jjl. 

The interface structure in the dimple region results 
from a balance between viscosity and surface tension. To 
describe the profile near the dimple position Zd, we there- 
fore look for a similarity solution of ([!]) 



h = C& 2/3 H (n) , ri= (z-z d )/Ca 1/3 . 



(3) 



where Q £ [0, Q max ]- The dimple with Q max corresponds 
to the lowest possible thickness (lower bound of grey re- 
gion in Fig. [2J, while the case Q — marks the upper 
boundary for hf with nonnegative flux. It is worth not- 
ing that Q — is only the upper bound in the case of the 
dip-coating geometry. If an injection of liquid is added 
at the top of the plate, one could match much thicker 
films to the bath. The structure of the films will become 
rather different when the downward flow is much larger 
than the plate velocity, for which the analysis crosses over 
to flow down a wall at rest [T3"] . 

The similarity solution of ([lj corresponding to the film 
region just above the dimple is 



Zd 



Ca 1 / 6 



(6) 



Following LLD, this solution has the additional property 



In (roll, the ft,-scale Ca 1 ^ 2 is dictated by the film thickness 
([H]) for small Ca. Thus the typical curvature h" oc Ca 1 / 6 
vanishes in the limit Ca <C 1, which establishes the 
boundary condition H'^ = quoted above. The slopes in 
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FIG. 3: (a) The three asymptotic regions of the dimple so- 
lutions, and the scaling of the respective film thickness with 
Ca. (b) The flux q as function of hf at a fixed Ca = 10 -4 , cf. 
(2|. The shaded region shows the range of thick solutions at 
that value of Ca. The horizontal line illustrates that the flux 
through the thick film is the same as through the thin dimple 
region. 




FIG. 4: Various solutions of Q: (a) Dimple profiles with Q = 
(solid), Q max = 1.376- • • (dashed), and the Landau-Levich 
film solution for which Q ~ 0.946 (dotted), (b) All curves 
have H'leo = v2 to match the curvature of the reservoir, and 
H'Je — to match to the film. 

the dimple and film regions, as given by ^ and (pj, re- 
spectively, both scale as Ca 1 ^ 3 and can thus be matched. 
To establish our central result ^ , a more detailed anal- 
ysis of the film solutions ([6| is not necessary. 

Note that the final approach onto the flat film involves 
(stationary) capillary waves on the film surface that are 
exponentially damped (see Fig. [3] (a)). From the point 
of view of the lubrication equation ([!]), these waves pro- 
vide an additional degree of freedom that is necessary to 
achieve different values of flux [U [21 [TS]. Linearizing 
around the flat film, one finds that such waves can in 
principle exist when hf > Ca 1 / 2 . As can be seen from 
([5]), however, this condition is not sufficient to explain 
the new films. Similar dimpled solutions have also been 
found for Marangoni-driven flows |16j . 

Experimental methods - In our experiment, a sili- 
con wafer is withdrawn vertically from a bath of sili- 
cone oil (viscosity r\ = 4.95 Pa-s, surface tension 7 = 



0.0203 N-m^ 1 , density p — 970kg-m~ 3 , molecular size 
70 nm) by a step-motor. The silicon wafer is totally 
wetted by the silicone oil. When coated with the fluo- 
rinated surfactant FC725, partial wetting conditions are 
obtained. Depending on the protocol, the contact angle 
lies between 48° and 55° , with an hysteresis of 5° . The 
film thickness hf is measured by spectrometry. A re- 
flection probe made of a tight bundle of six illumination 
fibers around one read fiber is placed at 5 mm from the 
plate. It is connected to a Tungsten Halogen light source 
and to a diffraction-grating spectrometer resolving visible 
and near-infrared wavelengths. The relative resolution, 
limited by the entrance slit and the size of the photo- 
sensitive elements, is around 0.25 nm. A 3648-element 
linear CCD-array detector measures the intensity as a 
function of the wavelength A, averaged over 4s. Defects 
of the CCD are calibrated in the absence of film and cor- 
rected to improve the signal-to-noise ratio. Once a film 
is present, the spectrum contains oscillations of the form 
cos(47rnft,y/A), where n — 1.4034 is the refractive index of 
silicone oil. By fitting the spectrum, the thickness hf can 
be measured to within 0.2% up to a thickness of 240 /im 
(cf. Fig. [5). 

The spatial profile of the entrained film is measured by 
placing a 200 /im wire at a distance d = 14 ± 1 mm from 
the silicon wafer. The mirror image of the wire (reflected 
in the silicon plate) is distorted by refraction through 
the liquid/ vapor interface of local slope h! , making the 
system equivalent to a wire placed at a distance d behind 
a prism of (small) angle 2h' . At small angles, the rays are 
thus deflected by an angle 2(n — l)h' in the direction of 
steepest slope, independent of the incident angle. Hence 
the mirror image of the wire is vertically displaced by a 
distance 2(n — l)dh' . The wire is imaged with a 2048 
x 2048 CCD camera, fitted with a 60 mm macro lens. 
The position of the wire image is determined by image 
inter-correlation, achieving sub-pixel resolution. 

We have calibrated the relation between the displace- 
ment and the local interface slope by two independent 
methods. First, we have replaced the oil film by op- 
tical glass prisms (n = 1.52) of angle 1°, 2°, 3° and 
4°. We obtain a resolution of 950 ± 30 pixels per unit 
slope (16.5 pixels/o), once the index difference between 
oil and glass is taken into account. Second, using the 
step-motor, we have determined an image resolution of 
11.7 /im/pixels, which leads to 970 ± 70 pixels per unit 
slope. To summarize, the displacement of the image of 
the wire gives the local interface slope from which we re- 
construct the film profile, by integration. The constant 
of integration is fixed in the flat region of the film, using 
the value of hf measured by spectrometry. There are two 
sources of error in the final reconstruction: an absolute 
uncertainty of 3% in the slope, and the error from the 
spatial integration of the random noise present for each 
point. 

Finally, we return to the question of how a particular 
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FIG. 5: Measurement of the film thickness by spectroscopy. 
The free charge carriers produced by photon absorption are 
accumulated during 0.4 s, counted and displayed as a function 
of the inverse wavelength A -1 . This spectrum is averaged 
over 100 realizations and corrected from the effect of CCD 
heterogeneities, calibrated previously. It is then split into the 
sum of a smooth spectrum, characteristic of the light source, 
and an oscillating spectrum, shown above (symbols). It is 
finally fitted by a sine modulated in amplitude by a Gaussian 
envelope (solid line), (a) LLD film at Ca = 9.75 x 10~ 2 . 
The best fit gives hf = 61.24 nm with (formally) a statistical 
accuracy of 0.005 nm. However, the reproducibility of two 
complete experimental runs is only to within 0.1 nm. (b) 
Thick film at Ca = 8.53 x 10~ 2 . The best fit gives hf = 
214.68 nm with (formally) a statistical accuracy of 0.05 nm 
and a reproducibility to within 1 nm. Thick films are at 
the limit of resolution, hence the interference pattern is only 
visible in the near-infrared range of wavelengths. 



type of film solution is realized experimentally. As shown 
in |8J, the front of a liquid film of partially wetting fluid 
possesses a characteristic speed of recession Ca/ relative 
to the substrate, which is essentially set by the contact 
angle. Since the flux q through the contact line is zero, 
according to ^ this translates into a characteristic film 
thickness hf = (30a/) 1 / 2 . Thus if the plate speed Ca 
is marginally above Ca/, the contact line moves up the 
plate very slowly, and the flux through the liquid film is 
close to zero: one has prepared the state (/i/,/i 2 /3) at 
the upper end of the thick film solutions. 

If the plate speed is now slowly increased, as we did 
in an experimental run shown in Fig. [2] by the arrows, 
one moves horizontally into the grey region, since hf is 
fixed by the motion of the contact line. As the end of the 
grey region is reached, the solution falls off the region of 
allowed solutions, and a LLD solution is realized instead 
(cf. Fig. |2j. During this transition, the dimple detaches 
from the bath, and moves up the plate, leaving behind 
the LLD film. In effect, this is the situation described 
in our earlier paper [5]. Note that we were not able to 



take a measurement right above the transition, since the 
LLD film takes too long to develop. Of course the above 
does not imply that thick films can only be produced in 
partial wetting situations. The most promising experi- 
mental setup to produce a thick film is that of viscous 
flow in the interior of a cylinder, which can be rotated in 
both senses. This offers considerably more flexibility in 
preparing thicker films. 

In conclusion, we have calculated a phase diagram for 
film coating from a bath. At a given speed, two differ- 
ent types of film solutions are possible. The novel, thick 
films described here are easily realized experimentally, 
and experiment agrees extremely well with theoretical 
prediction. 
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